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Abstract 

It is shown that the temperature dependence of Langmuir constants contains all the informa- 
tion needed to determine spherically averaged intermolecular potentials. An analytical "inversion" 
method based on the standard statistical model of van der Waals and Platteeuw is presented which 
extracts cell potentials directly from experimental data. The method is applied to ethane and cy- 
clopropane clathrate-hydrates, and the resulting potentials arc much simpler and more meaningful 
than those obtained by the usual method of numerical fitting with Kihara potentials. 

1 Introduction 

A central mission of chemical physics is to determine intermolecular interactions from experimental 
phase equilibrium data. This is generahy a very difficult task because macroscopic equilibrium 
constants reflect averaging over vast numbers of poorly understood microscopic degrees of freedom. 
Since the advent of the computer, empirically guided numerical fitting has become the standard 
method to obtain intermolecular potentials. An appealing and often overlooked alternative, however, 
is to solve inverse problems based on simple statistical mechanical models. 

A well-known class of inverse problems is to determine the density of states n{E) from the 
partition function Q{(3), or one of its derivatives, as a function of inverse temperature /? = 1/kT. For 
classical Maxwell-Boltzmann statistics, these quantities are simply related by an Laplace transform |j|] 

Q{(3) = / e-^'n{e)de (1) 

^0 

which is straightforward to invert. Analogous quantum- mechanical problems have also been solved 
using Laplace (or Mellin) transform methods. For example, in Bose-Einstein statistics, 

the phonon density of states can be obtained from the specific heat of a crystal ||2|, ^, ^, ^, ^, ^ ^, 
and the area distribution of a blackbody radiator can be obtained from its power spectrum |jlO|, 
12, ^, Similarly, in Fermi-Dirac statistics, 



1 



the band structure of a doped semiconductor can be obtained from the temperature dependence of 
its carrier density ||l^. In each of these examples, it is possible to extract the microscopic density 
of states from a temperature-dependent thermodynamic quantity because each is a function of just 
one variable. 

In the more complicated case of chemical systems, thermodynamic quantities are related to 
classical configurational integrals, 

-/3*{r), 



Z{P) = / e-^^^'^^df, (4) 
Jv 

where the total energy e in Eq. (0) is replaced by the intermolecular potential <I>(r) , and the 
integral extends over the interaction volume V. Unfortunately, these integrals are multi-dimensional, 
which leads to under-determined inverse problems. Perhaps as a result, statistical inversion methods 
have apparently not been developed for intermolecular potentials in chemical systems, and instead 
empirical fitting has been used exclusively to describe phase equilibrium data. 

In this article, we show that within the common spherical-cell approximation the intermolecular 
potential is completely determined by the temperature dependence of the Langmuir constant. In this 
case, the linear integral equation (||) is replaced by a non-linear equation of the form 

roo 

Z{I3) = 47r / e-'^"'(''V2dr (5) 
where w{r) is a spherically averaged "cell potential" |24]. As shown below in the important case of 



Langmuir constants for clathr ate- hydrates, this simple inverse problem can be solved exactly without 
resorting to numerical fitting schemes. 

Before proceeding, we mention some related ideas in the recent literature of solid state physics. 
As with chemical systems, empirical fitting is also the standard approach to derive interatomic poten- 
tials for metals and semiconductors. Since the pioneering work of Carlsson, Gelatt, and Ehrenreich 



in 1980 [14, 15 1, however, exact inversion methods have been developed to obtain potentials from co- 



hesive energy curves ||5|, |T^, 17, 18, 19, |2^. These theoretical advances, discussed briefly in Appendix 
A, have recently led to improvements in the modeling of silicon, beyond what had previously been 
obtained by empirical fitting alone [20, |2l[. Inspired by such developments in solid state physics, 



here we seek similar insights into clathrate-hydrate intermolecular forces, albeit using a very different 
statistical mechanical formalism. 



2 The Inverse Problem for Clathrate-Hydrates 
2.1 The Statistical Theory of van der Waals and Platteeuw 

Clathrate-hydrates exist throughout nature and are potentially very useful technological materi- 
als |2^. For example, existing methane hydrates are believed to hold much more energy than any 
fossil fuel in use today. Carbon dioxide hydrates are being considered as effective materials for the 
sequestration and/or storage of CO2. In spite of their great importance, however, the theory of 
clathrate-hydrate phase behavior is not very well developed, still relying for the most part on the 
ad hoc empirical fitting of experimental data. Therefore, we have chosen to develop our statistical 
inversion method in the specific context of clathrate-hydrate chemistry. 

Since being introduced in 1959, the statistical thermodynamical model of van der Waals and 
Platteeuw (vdWP) has been used almost exclusively to model the phase behavior of clathrate- 
hydrates, usually together with a spherical cell (SC) model for the interaction potential between 
the enclathrated or "guest" molecule and the cage of the clathrate-hydrate. The SC model was 
also introduced by vdWP, inspired by an analogous approximation made by Lennard-Jones and 



Devonshire in the case of liquids [^, 24 1 . 
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In the general formulation of vdWP |23|, the chemical potential difference between an empty, 
unstable hydrate structure with no guest molecules, labeled MT, and the stable hydrate, labeled H, 
is related to the so-called Langmuir hydrate constant Cji and the fugacity of the guest molecule fj 

Aa.^^^-^ = kTY,vdr^{l + Y, Cjjj) (6) 

i J 

where i designates the type of cage, Ui the number of cages of type i per water molecule and J 
the type of guest molecule. In practice, experimental phase equilibria data is used to determine 

The connection with intermolecular forces within vdWP theory is made by expressing the Lang- 
muir hydrate constant as the configurational integral Zji divided by kT, which is written explicitly 
as an integral over the volume V 

CjAT) = I e-^^~'''^^^''''^'^^/^'^r'^sinesinidrded<t)dadid-i (7) 

S'K'^kT Jy 

where $(r, 0, (/>, a, ^, 7) is the general six dimensional form of the interaction potential between the 
guest molecule at spherical coordinates (r, 6, cj)) oriented with Euler angles (a, ^, 7) with respect to 
all of the water molecules in the clathrate-hydrate. 

In the SC approximation, which is made without any careful mathematical justification, the 
intermolecular potential <I> is replaced by a spherically averaged cell potential w{r), which reduces 
the Langmuir constant formula (0) to a single, radial integration 

CMT) = ^j\-^^y'^r'dr (8) 

where the cutoff distance R is often arbitrarily taken as the radius of the cage. (As shown below, the 
exact value rarely matters because temperatures are typically so low that the high energy portion 
of the cage r ~ i? makes a negligible contribution to the integral.) Although the SC approximation 
may appear to be a drastic simplification, it is nevertheless very useful for theoretical studies of 
intermolecular forces based on Langmuir constant measurements. 

2.2 Numerical Fitting Schemes 

Before this work, the functional form of the cell potential w{r) has always been obtained by first 
choosing a model interaction potential between the guest molecule in a cage and each nearest neigh- 
bor water molecule essentially ad hoc, and then performing the spherical average from (^) to (^) 
analytically. The most common potential form in use today is the Kihara potential, which is simply 
a shifted Lennard-Jones potential with a hard-core. Using the Kihara potential and spherically av- 
eraging the interaction energy, typically over the first-shell only, yields the following functional form 
for w{r): 

-(0 = 2-[|^?;:(^io + 1^11) - ^(^4 + ^s,)] (9) 

where 

and z is the coordination number, R is the radius of the cage, and a, e, and a are the Kihara 
parameters. As a result of the averaging process leading from (Q) to @, the functional form of w{r) 
is fairly complicated, and the parameters e and a are generally determined by fitting monovarient 



equilibrium temperature-pressure data numerically [P2| , 25] 



3 



There are several serious drawbacks to this ubiquitous numerical fitting procedure, which suggest 
that the Kihara parameters lack any physical significance: (i) The Kihara parameters are not unique, 
and many different sets can fit the experimental data well; (ii) the Kihara parameters found by fitting 
Langmuir curves do not match those of found by fitting other experimental data, such as the second 



virial coefficient or the gas viscosity [22|; and {iii) comparisons of Langmuir constants found via 
the SC approximation ^ and via explicit multi-dimensional quadrature (0) show that the two can 
differ by over 12 orders of magnitude ^] (which results from the exponentially strong sensitivity 
of the Langmuir constant to changes in the cell potential). These problems call into question the 
validity of using the Kihara potential as the basis for the empirical fitting, and even the use of the 
SC approximation itself. 

2.3 Inversion of Langmuir Curves 

It would clearly be preferable to extract more reliable information about the interatomic forces 
in clathrate-hydrates directly from experimental data without any ad hoc assumptions about their 
functional form. In principle, such an approach is possible for clathrate-hydrates which contain a 
single type of guest molecule occupying only one type of cage. In this case, each of the sums in 
Eq. (^) contains only one term, and by using an equation of state to compute the fugacity /j, the 
Langmuir constant Cji can be determined directly from experimental phase equilibria data. Typical 
data sets obtained in this manner are shown in Fig. |l] for Structure I ethane and cyclopropane 



clathrate-hydrates [28| 



Because the full potential ^ in (^) is multi-dimensional (while the Langmuir constant only de- 
pends on a single parameter T), the general vdWP theory is too complex to pose a well-defined 
inverse problem for the interatomic forces. The SC approximation, on the other hand, introduces a 
very convenient theoretical construct, the spherically averaged potential w{r), which has the same 
dimensionality as the Langmuir curve Cji(T) of a single type of guest molecule occupying a single 
type of cage. (Since we consider only this case, we drop the subscripts Ji hereafter.) Although one 
can question the accuracy of the SC approximation, its simplicity at least allows precise connections 
to be made between the Langmuir curve and the cell potential. 

As an appealing alternative to empirical fitting, in this article we view Eq. (^) as an integral 
equation to be solved analytically for w{r), given a particular Langmuir curve C{T). Letting (3 = 
1/kT, we rewrite (|8|) as 

roo 

C{(3) = 47r/3 / e-'3'"(''V2dr, (11) 





where we have also set the upper limit of integration to R = oo, which introduces negligible errors due 
to the very low temperatures (large (5) accessible in experiments. (This will be justified a posteriori 
with a precise definition of "low" temperatures below.) Note that Eq. ( pT|) has the form of Eq. (^) 
with Z{P) = pC{P). 

2.4 Application to Experimental Data 

In our analytical approach, some straight-forward fitting of the raw experimental data is needed to 
construct the function C(/3), but after that, the "inversion" process leading to w{r) is exact. For 
example, typical sets of experimental data are well described by a van't Hoff temperature dependence 

Cif3) = Coe^l" (12) 

as shown in Fig. ^ for ethane and cyclopropane clathrate hydrates [^], and the constant m is 
generally positive. (Note that the exponential dependence which we call "van't Hoff dependence" in 
this paper can be expected based on quite general thermodynamic considerations [p9|].) 
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Aided by the analysis, the quahty and functional form of these fits are discussed below in section |^. 
In order to allow for deviations from the dominant van't Hoff behavior, however, in this article we 
consider the more general form 

C{P) = pF{P)e"'f^ (13) 

where m is a constant defined by 

m = lim logC{p)/p (14) 

/3^oo 

whenever this limit exists and is finite, i.e. when the prefactor F{f3) in Eq. ( |T^ ) is dominated by the 
exponential term at low temperatures. We exclude the possibility of hyper-exponential behavior at 
low temperatures, logC//3 oo, which is not physically meaningful, as explained below. The set of 
possible prefactors includes power-laws, F{/3) = f5~^, as well as various rational functions. 

The rest of article is organized as follows. In section ^ we discuss various necessary and sufficient 
conditions for the existence of physically reasonable solutions, and we also derive the asymptotics of 
the Langmuir curve at low temperature from the behavior of the cell potential at near its minimum. 
In section we perform the analysis in the general case, and in section ^, we discuss the specific 
case of van't Hoff dependence ([l2|), which leads to a cubic solution as well as various unphysical 
solutions involving cusps. In section |6|, we derive analytical solutions for several different temperature 
dependences, which reveal the significance of possible deviations from van't Hoff behavior for the form 
of the potential w{r). The theoretical curves are compared with the experimental data in section |^, 
and the physical conclusions of the analysis are summarized in section ^. Relevant mathematical 
theorems are proved in the Appendix B. 



3 General Analysis of the Inverse Problem 
3.1 Necessary Conditions for the Existence of Solutions 

On physical grounds, it expected that the cell potential w{r) is continuous (at least piecewise) and 
has a finite minimum at ro > somewhere inside the clathrate cage, w{r) > w{ro) = Wq- We also 
allow the possibility that w{r) is infinite for certain values of r (e.g. outside a "hard- wall radius") 
by simply omitted such values from the integration in Eq. ([ll|). As proved in the Appendix B, these 
simple physical requirements suffice to imply the asymptotic relation (14), where m = —Wq- They 



also place important constraints on the prefactor F{0) defined in (^), which must be 

(i) analytic in the half plane Re/3 > c and 

(ii) real, positive and non-increasing for /? > c on the real axis 

where c > is a real number. (Note that we view the inverse temperature /? as a complex variable, 
for reasons soon to become clear.) Moreover, if the set 

S = {r> Q\wo < w{r) < oo} (15) 

has nonzero measure, then F{P) is strictly decreasing on the positive real axis. It is straightforward to 
generalize these rigorous results to the multi-dimensional integral of vdWP theory, Eq. (|^, without 
making the spherical cell approximation, as described in the Appendix B, but hereafter we discuss 
only the spherically averaged integral equation, Eq. ([TTl) , because it makes possible an exact inversion. 

In this section, we give simple arguments to explain the results proved in Theorem 1 of the 
Appendix B. First, we consider the illustrative example of a constant cell potential with a hard wall 
at r = Vhw > 0, 

, . i Wo if < r < rhw /-i^n 
' ' oo if r > r/,^ 
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which satisfies the assumptions stated above. The integral (11) is easily performed in this case to 
yield 

C(/3) = ^vrrL/3e-°^ (17) 

which implies m = —Wo, the well depth, and F{P) = ^Trrf^^, the volume of negative energy. Con- 
sistent with the general results above, F{P) is constant in this case since S" = 0. For continuous 
potentials w{r), however, the prefactor F{P) must be strictly decreasing because 5 7^ 0. 

The integral equation (|l^) can be simplified by a change of variables from radius to volume. In 
terms of a shifted cell potential versus volume 

47r 

u{x) = w{r) — Wo where x = — r (18) 

3 

the integral equation is reduced to the form 

C(/3) = (3F{l3)e-'"^'''^^ (19) 

where 

roo 

F{P) = / e-^"(^)dx. (20) 
Jo 

Since u{x) > by construction, the function F{P) is clearly non-increasing. In the Appendix B, it 
is proved that if the potential varies continuously near its minimum (in a very general sense), then 
F{P) does not decay exponentially. Since F{f3) also positive and bounded above, we conclude 

lim logC(/3)//3 = -w;o (21) 

which implies m = —Wq- Therefore, the slope of the Langmuir curve on a "van't Hoff plot" (logC 
versus /3 = 1/kT) in the low temperature limit is equal to (minus) the minimum energy of the cell 
potential. Since it is generally observed that m is positive, as in the case of ethane and cyclopropane 
clathrate-hydrates shown in Figs. ^|-|2|, the cell potential must be attractive, Wo = —m < 0, which 
simply indicates that the total internal energy is lowered by the introduction of guest molecules into 
clathrate-hydrates. 

The fact that F{P) must be nondecreasing has important consequences for the existence of 
solutions which are piecewise continuous and bounded below. For example, consider the class of 
Langmuir curves of the form 

C{P) = p^e""^ (22) 

which is useful in fitting experimental data (see below). We have already addressed the borderline 
case, = 1, in which a discontinuous hard- wall solution is possible. Although it is not obvious 
a priori, there are no solutions to the inverse problem if > 1, since in that case F{/3) = /3'^^^ is 
increasing. On the other hand, if < 1, then well-behaved continuous solutions are possible, because 
F{f3) is strictly decreasing. 



3.2 Low Temperature Asymptotics of the Langmuir Curve 

/,Prom Eq. (^l]), we see that the minimum energy Wq = w{ro) determines the leading order asymp- 
totics of C(/3) in the low temperature limit. More generally, one would expect that C(/3) at low 
temperatures is completely determined by shape of the cell potential at low energies, close to its 
minimum. Using standard methods for the asymptotic expansion of Laplace integrals [^], it is 
straightforward to provide a mathematical basis for this intuition. For simplicity, here we consider 
the usual case of a parabolic minimum 

w{r) =Wo + h{r - rof + 0((r - rof) (23) 
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for some constants k > and > 0, although below we will derive many exact solutions with 
non-parabolic minima. 

Due to the factor of appearing in the integrand in Eq. (11), the two cases of a non-central 



or central minimum, > and Tq = 0, respectively, must be treated separately. Physically, this 
qualitative difference between central and non-central-wells is due to the spherical averaging process 
going from Eq. (0) to Eq. ([Tl|): A central well in w{r) corresponds to a unique local minimum of the 
multidimensional potential <I>, but a non-central-well in w{r) corresponds to a nonlocal minimum of 
$ which is smeared across a sphere of radius Tq. 

Beginning with non-central-well case, r > 0, we have the following asymptotics as Re/? — oo: 

C{/3) ~ AnrlP p^\-Pi''"o+\k{r-r,f)^^ 



47rr^ 



{^f-'fe-^"' (24) 



O I 

which is the usual leading order term in the expansion of a Laplace integral |3^. Therefore, the ex- 
perimental signature of a non-central- well is a Langmuir constant which behaves at low temperatures 
like 

C{p) ~ Coe™^/3^/2 as oo. (25) 



Comparing (24) and (P5D, we can identify the well depth Wq = —m, consistent with the general 
arguments above, but it is impossible to determine independently the location of ro and the curvature 
k of the minimum. Instead, any and k satisfying Airr'^^/^jr/Tc = Co would exactly reproduce the 
same large-/? asymptotics of the Langmuir curve (as would a completely different central- well solution 
described in section |6|). This degeneracy of non-central- well solutions revealed in the low temperature 
asymptotics is actually characteristic of all non-central-well solutions, as explained below. 

In the central-well case, = 0, the asymptotics must be carried out more carefully because the 
leading term derived in (|2^) vanishes: 



C7(/3) ~ 47r/3 / g-^^^'^+^'^^'V^dr 

1/2 







/ 9 \ i/^ /-oo 



The experimental signature of a parabolic central well in the Langmuir curve, 

C(/3) ~ ae"^/3-^/2 as /? ^ oo, (27) 



is qualitatively different from (pSD, which provides an unambiguous way to separate the two cases 
using low temperature measurements. Moreover, unlike the non-central-well case, the curvature k = 

—2/3 

27rCo of a parabolic central minimum is uniquely determined by the low temperature asymptotics 
of the Langmuir curve. Consistent with asymptotic results, we shall see in section ^ that central-well 
solutions to the inverse problem are unique, while nan- central- well solutions are not. 

3.3 Sufficient Conditions for the Existence of Solutions 

The primary difficulty in solving Eq. (0) lies in its being a nonlinear integral equation of the "first 



kind" for which no general theory of the existence and uniqueness of solutions exists [31, 32]. In the 
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linear case, however, there is a special class of first-kind equations which can be solved using Laplace, 
Fourier or Mellin transforms, namely integral equations of the additive or multiplicative convolution 
type H 



^{x) = / Kix-y)i^{y)dy (28) 

^— oo 

or ^{x) = / K{xy)i:{y)dy, (29) 



respectively, where 'iIj{x) is the unknown function and ^{x) is given. Integral equations of the 
multiplicative form ( p9D often arise in statistical mechanics as explained in the Introduction. 

Although our nonlinear, first-kind equation (|l^) is not of the convolution type because the un- 
known function w{r) appears in the exponent, it does somewhat resemble a Laplace transform. This 
connection is more obvious in the alternative formulation (|20|) relating -F(/3) and u{x), which is 



equivalent to the original equation (11) according to the analysis above. In the next section, it is 
shown that physically reasonable solutions exist whenever F{j3) has an inverse Laplace transform 
f{y) which is positive, nondecreasing and non-constant for y > 0. (Sufficient conditions on F{I3) to 
ensure these properties of f{y) are given in Appendix B.) 



4 Analytical Solutions for Arbitrary Langmuir Curves 
4.1 The Unique Central- Well Solution 

It is tempting to change variables y = u{x) in the integral (^) to reduce it to a Laplace transform, 
but care must be taken since u^^(y) may not be single- valued. This leads us to treat solutions 
which are monotonic separately from those from those which are not, an important distinction fore- 
shadowed by the asymptotic analysis above. As a natural first case, we seek differentiable solutions 
u{x) which are strictly increasing without bound (u(oo) = oo) from a central minimum {u{0) = 0). 
Such "central- well solutions" correspond to cell potentials w{r) which are strictly increasing from 
a finite minimum w{0) = Wq at the center of the cage. We proceed by considering the inverse cell 
potential v{y) = u~^{y) with units of volume as a function of energy, which is single- valued and 
strictly increasing with f (0) = 0, as shown in Fig. 0(a). With the substitution y = u{x), Eq. (pO|) is 



= / -^r-T7^ = / e-''yv'iy)dy. (30) 



reduced to Laplace's integral equation |35| for the unknown function v'{y) 

coo c-Pvd^i roo 
lo u'{u-^{y)) Jo 

Upon taking inverse Laplace transforms, we arrive at a differential equation for v{y) 

v'{y) = f{y) (31) 

whose unique solution is 

v{y)=rf{y)dy (32) 







using the boundary condition f (0) = 0. According to (pjl|), the continuity of f{y) for y > (which 
is not assumed) would the guarantee differentiability of v{y) for y > 0, and hence of u{x) for x > 0. 
Equivalently, we can also simplify (|30|) with an integration by parts 



oo 



F{[j) = P / e-"yv{y)dy. (33) 
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cm = ^ = (35) 



Therefore, the inverse ceh potential is given by 

viy) = 9{y) (34) 
where g{y) is the inverse Laplace transform of the function 

F(/3) _ C(/3)e ^"'° 

The cell potential u{x) is determined implicitly by the algebraic equation 

g{u) = X. (36) 
Returning to the original variables, we have a general expression for w{r) in the central-well case 

w{r) = Wo + g-^ (^-Kr^^ . (37) 

This equation uniquely determines the central-well potential that exactly reproduces any admissible 
Langmuir curve. 

4.2 Non-Central- Well Solutions 

The simplest kind of non-central-well solution is the central-well (^^ shifted by a "hard-core" radius 
The > 

{oo if < r < 

nj, + g-^[ln{r^-rl)] if r > r,, ^^^^ 

which exemplifies a peculiar general property of our integral equation: An arbitrary hard core can 
be added to any solution. Note that, if u{x) is any solution of the rescaled equation (|20|), then so is 



u{x) = l , if0^x<x,^ 

\ U{X- Xhc) X>Xhc 



for any hard-core volume Xhc ^ 0. The proof is simple: 



oo roc 



^-(3u(x)d^= g-/3n(x-x„,)^^ ^ / e-'^"(^)dx = F(/3). (40) 

J Xhc Jo 

Physically, a hard core for the cell potential could represent the presence of a second guest molecule (in 
a spherically symmetric model) in the same clathrate-hydrate cage. Alternatively, a hard-core could 
represent a water molecule (again in a spherically symmetric model) at the node of several adjacent 
clathrate cages, in which case the cell potential actually describes the "super-cage" surrounding the 
central water molecule. 

The arbitrary hard-core just described only hints at the vast multiplicity of non-monotonic so- 




lutions to the integral equation (20), which is a common characteristic of first-kind equations 
Next, we consider the general case of a non-central- well, shown in Fig. |3|(b), which includes (^8]) as 
a special case. To be precise, we seek continuous solutions u{x) on an interval xi < x < X2 com- 
posed of a non-increasing function u-{x) and a nondecreasing function u-{x) which are piecewise 
differentiable and non-negative. We also allow for a possible hard-core in the central region x < xi 
as well as a "hard wall" beyond the clathrate cage boundary x > X2- The general form of such a 
non-central-well solution is 

if < X < xi 

if X < X2 
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where U-{xo) = u+{xo)- We do not assume u'_^{xo) = u'_{xo), which would imply differentiability 
at the minimum u'{xo) = 0, although we do not rule out this case either. Instead, we allow for the 
mathematical possibility of a discontinuous first derivative at Xo, i-e. a "cusp" at the well position, 
at least for the moment. 

As before, it is convenient to express the solution (|4l|) in terms of two differentiable functions 
V-{y) = uZ^{y) and v^{y) = u'^^{y) which describe the multi- valued inverse cell potential. Note 
that f-_(oo) = xi, v^{oo) = X2 and f-(0) = f+(0) = Xq- In terms of the inverse cell potentials, the 
integral equation ( pO|) takes the form 

F{(j) = I e-^"(^)dx 



which implies 







oo 

e 





v'^{y)-v'__{y) = f{y). (43) 

In this case, the continuity of /(y) would only guarantee the differentiability of the difference f+(y) — 
V-{y), but not of the individual functions f+(?/) and V-{y). Integrating (|4^ ) by parts before taking 
the inverse transform yields a general expression for the solution 

v+{y) - v^{y) = g{y) (44) 

where again g{y) is the inverse Laplace transform of F{(3)/(3. Unfortunately, we have two unknown 
functions and only one equation, so the set of non-central-well solutions is infinite. 

The scaled Langmuir curve F{(3) uniquely determines only u+(y) —V-{y), the volume difference as 
a function of energy between the two branches of the scaled cell potential u{x), but not the branches 
V-^{y) and V-{y) themselves. An infinite variety of non-central- well solutions, which exactly reproduce 
the same Langmuir curve as the central well solution, can be easily generated by choosing any non- 
increasing, non-negative, piecewise differentiable function V-{y) such that the function w+(y) defined 
by ( ^4| ) is nondecreasing. Even the position of the well V- (0) = Xo can be chosen arbitrarily. 

For example, one such family of solutions with a central "soft-core" (xi = 0) is given by 

if0^x<a;. 

it X > 




where 



and 

v+{y) = ul\y) = v^{y) + g{y), (47) 

for any a,b > and Xo > 0. (In the limit a — > 0, we recover the unique central-well solution.) Note 
that yc = m(0) = {xo/aY^^ is the height of the central maximum of the potential. These solutions, 
all derived from a single Langmuir curve, exist whenever g{y) increases quickly enough that f+(y) is 
nondecreasing, which is guaranteed if g'{y) > ahy^"^ for < y < yc- 

Another family of non-central-well solutions with a soft-core can be constructed with the choice 

I if y > yc ^ ^ 
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for any < A < 1 and Xo > 0, where Vc = 9 ^(a^o/A). In this case, the cell potential is easily 
expressed in terms of g~^{x) as 

/ g-i ^xo^^ if < X < Xo 
u{x) = I g-^ (f^) if Xo < X < Xe (49) 
[ g^^{x) if X > Xc 

where Xc = v+{yc) = Xo/A. This class of solutions exists whenever g{y) is nondecreasing (or f{y) > 0). 
If yc = oo, then there is a hard core m(0) = oo. Otherwise, if there is a soft core u(0) = yc < oo, 
then there is typically a cusp (discontinuous derivative) at Xc, as explained below. 

As demonstrated by the preceding examples, it is simple to generate an enormous variety of 
non-central-well solutions, with an arbitrarily shaped soft or hard core, and an arbitrary position of 
the minimum. In spite of the multiplicity of non-central-well solutions, however, our analysis of the 
inverse problem at least determines f+(?/) — V-{y) uniquely from any experimental Langmuir curve. 
This important analytical constraint is not satisfied by empirical fitting procedures. 

4.3 Soft Cores and Outer Cusps 

Non-central- well solutions with a soft-core satisfy f _ (y) = for y > yc > 0, as in the examples above. 
In such cases, i'+(y) = g{y) for y > yc regardless of whether or not there is an outer hard wall, which 
implies that u(x) = g~^{x) for x > Xc, where Xc = g~^{yc)- If fiu) is continuous for y > 0, then, 
unless U-{x) has an "inverted cusp" at the origin {v'_{yc) = and u'_(0^) = — oo), any non-central- 
well solution u{x) with a soft-core must have a cusp at x = Xc, as in the examples above. This "outer 
cusp" in u{x) could only be avoided if /(y) itself has a cusp at yc which would allow f+(y) to be 
continuous. However, an inverted cusp in u{x) at the origin does not necessarily imply an cusp in 
w{r) at the origin due the transformation x = l^rr^. For example, if V-{y) ~ (yc — y)^^^ as y ^ yc, 
or n(x) ^ yc — x"^^^ as x ^ 0, then w{r) would have a physically reasonable, parabolic soft core 
w{r) ^ Wo + yc — (47r/3)^''^r^ as r ^ 0. Nevertheless, even in such cases, if /(y) were continuous for 
all y > 0, then both u{x) and w{r) would have unphysical second-derivative discontinuities at x = Xc 
related to the soft core. In general, a continuously differentiable, non-central-well solution with a 
central soft core could only arise if /(y) were discontinuous at some yc > 0, and such discontinuities 
are generally not present. 

4.4 Cusps at a Non-Central Minimum 

As mentioned above, the behavior of the cell potential near its minimum (whether central or not) 
is determined by the behavior of the Langmuir curve at low temperature, or equivalently, at large 
inverse temperature, (3 = 1/T. The Laplace transform formalism makes this connection transparent 
and mathematically rigorous. The asymptotic behavior of G{f3) as Re/3 — > oo is related to the 
asymptotics of the inverse transform g(y) as y ^ 0, which in turn governs the local shape of the 
energy minimum through Eq. (^) for a central well or Eq. (|4^) for a non-central well. The leading 
order asymptotics has already been computed above for parabolic minima, but the general solutions 
above show how various non-local properties of the potential are related to finite temperature features 
of the Langmuir curve. Here, we comment on a subtle difference in differentiability between central 
and non-central- wells, related to the small-y behavior of /(y). 

For typical sets of experimental data, including the van't Hoff form (|l^), the prefactor F(P) has 
a bounded inverse Laplace transform in the neighborhood of the origin 

lim / (y) = /(O) < oo. (50) 
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This generally implies the existence of a cusp at a non-central minimum of u{x), which is signified 
by a nonzero right and/or left derivative. When u{x) is differentiable at its minimum, it satisfies 
u'_{x~) = n'_,_(x+) = 0. In the central-well case x^ = 0, the existence of a cusp in u{x) follows from 

m 

n'(0+) = l/v'{Q+) = l//(0) > 0, (51) 

but this does not imply a cusp in the unsealed potential w{r) as long as /(O) > because in that 
case 

w{r) -Wo = n(47rrV3) ~ (47r/3/(0))r3 as r ^ 0+. (52) 
In the non-central-well case, however, the bounded inverse transform (|50|) implies a cusp at the 



minimum because, v'_^_(0) — v'_{0) = /(O) < oo from (44) along with ^^1.(0) > and v'_(0) < implies 
that w+(0) < oo and/or v'_{0) > — oo which in turn implies u'^{0) > and/or u'_{0) < 0. Unlike the 
central-well case, however, a cusp in u{x) at the non-central minimum Xq > implies a cusp at the 
corresponding non-central minimum of w{r). Therefore, we conclude that whenever (50) holds, the 
only physically reasonable solution is the central- well solution (37). 



4.5 Asymptotics at High Energy and Temperature 

The high energy behavior of the cell potential is related to (but not completely determined by) the 
high temperature asymptotics of the Langmuir hydrate constant, through the function g{y)- For 
example, the cell potential would have a hard wall at X2 < oo, if and only if g{y) were unbounded 

lim g(y) = oo. (53) 

Since the empirical modeling of Langmuir curves using Kihara potentials assumes an outer hard- 
core, Eq. ( |53| ) could be used to test the suitability of using the Kihara potential form, although 
experimental data is often not available at sufficiently high temperatures to make a fully adequate 
comparison (see below). Whenever ( |53[ ) holds, the non-central- well solutions u{x) are also universally 
asymptotic to the central-well solution 

u{x) ~ g'^{x) (54) 

at large volumes x oo. This follows from (^^ and the fact that V-{y) is bounded, which implies 
V-i-{y) ~ g{y)- The exact inversions performed in section ^ provide further insight into the relationship 
between small /? asymptotics of the Langmuir constant and high energy behavior of the cell potential. 



5 Langmuir Curves with van't HofF Temperature Dependence 

Experimental Langmuir hydrate-constant curves C(/3) are well fit by an ideal van't Hoff temperature 
dependence (|l2|), demonstrated by straight lines on Arrhenius log-linear plots 

logC7 = m/3 + logC7o (55) 

as shown in Figs. |l| and ^ for ethane (Co = 4.733 x 10^^ atm^^, m = 9.4236 kcal/mol) and cy- 
clopropane (Co = 1.9041 X 10^'^ atm^^, m = 10.5939 kcal/mol) clathrate-hydrates [^8|. This data 
is analyzed carefully in section ^, where alternative functional forms are considered. In the ideal 
van't Hoff case, we have F{P) = Cq/P and G(/3) = Co/(3'^- The inverse Laplace transforms of these 
functions are simply f{y) = CoH{y) and g{y) = CoyH{y), respectively, where H{y) is the Heaviside 
step function. 
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We begin by discussing the unique central-well solution, which is illustrated by the solid line in 
Fig. ^ for the case of ethane. The central- well solution is linear in volume u{x) = g{x) = CoyH{y), 
and cubic in radius 

w{r) = - m. (56) 

A curious feature of this exact solution is that it has a vanishing "elastic constant", w"{0) = 0, a 
somewhat unphysical property which we address again in section |^. 

The simple form of (^) makes it very appealing as a means of interpreting experimental data 
with van't Hoff temperature dependence. We have already noted that the slope of a van't Hoff (Fig. 
2) plot of the Langmuir constant is equal to the well depth m = —Wq, but now we see that the 
y-intercept log Co is related to the well-size, e.g. measured by the volume of negative energy mCo- 
This volume corresponds to a spherical radius of 

which is 0.4180 A for ethane and 0.3208 A for cyclopropane. Since the van der Waals radius of ethane 
is less than that of cyclopropane, it makes physical sense that rf'^'^^^ > T-cydopropropane _ Moreover, 
these volumes fall within the ranges determined from two different experimental modeling approaches: 
Using the radius of the water cage from x-ray scattering experiments ||2^, Lennard- Jones potentials 
from gas viscosity data give 0.79 A for ethane and 0.61 A for cyclopropane ||3^, while computations 
with van der Waals radii give 0.18 A for ethane and 0.03 A for cyclopropane [p2| . 

There are infinitely many non-central-well solutions reproducing van't Hoff temperature depen- 
dence, but each of them has unphysical cusps (discontinuous derivatives). There will always be a 
cusp at the minimum of the potential, since f{y) satisfies the general condition (|50|). For example, 
the central-well solution can be shifted by an arbitrary hard-core radius > 



oo if < r < To 

"soT 



U!{r) = { 47r(r3-rg) . (58) 

— i oz — ui II r > To 



In the case of a soft core, there must be a second cusp in the outer branch of the potential at the 
same energy as the inner core due to the continuity of /(y), as explained above. This is illustrated 
by the following piecewise cubic family of soft-core solutions of the general form (E^) : 



87r|rg-r3| 



which are shown in Fig. Q in the case of ethane guest molecules. An infinite variety of other piecewise 
differentiable solutions exactly reproducing van't Hoff dependence of the Langmuir curve could easily 
be generated, as described above, but each would have unphysical cusps. 

Previous studies involving ad hoc fitting of Kihara potentials have reported non-central- wells [p^, 
but these empirical fits may be only approximating various exact, cusp-like, non-central-well solu- 
tions, such as those described above. Moreover, given that the central-well solution ( ^6[ ) can perfectly 
reproduce the experimental data, it is clear that the results obtained by fitting Kihara potentials 
to Langmuir curves are simply artifacts of the ad hoc functional form, without any physical signifi- 
cance. Kihara fits also assume a hard wall at the boundary of the clathrate cage (by construction), 
whereas all of the exact analytical solutions (both central and non-central-wells) have the asymptotic 
dependence 

471-7-3 

w{r) ~ ^ (60) 



if < r < 
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as r — i- oo according to (^). Any deviation from the cubic shape at large radii, such as a hard wall, 
would be indicated by a deviation from van't Hoff behavior at high temperatures, but such data 
would be difficult to attain in experiments (see below). 

The preceding analysis shows that the only physical information contained in a Langmuir curve 
with van't Hoff temperature dependence is the depth Wq and the effective radius r<j of the spherically 
averaged cell potential, which takes the unique form (|5^) in the central-well case. In hindsight, 
the simple two-parameter form of the potential is not surprising since a van't Hoff dependence is 
described by only two parameters, m and Co- It is clearly inappropriate to fit more complicated ad 
hoc functional forms, such as Eq. (P) derived from the Kihara potential, since they contain extraneous 
fitting parameters and do not reproduce the precise shape of any exact solution. 

6 Analysis of Possible Deviations from van't Hoff Behavior 
6.1 Dimensionless Formulation 

The general analysis above makes it possible to predict analytically the significance of possible 
deviations from van't Hoff temperature dependence, which could be present in the experimental 
data (see below). We have already discussed the experimental signatures of various low and high 
energy features of the cell potential in the asymptotics of the Langmuir curve. In this section, we 
derive exact solutions for Langmuir curves of the form (^) where F{f5) is a rational function. Such 
cases correspond to logarithmic corrections of linear behavior on a van't Hoff plot of the Langmuir 
curve, which are small enough over the accessible temperature range to be of experimental relevance, 
in spite of the dominant van't Hoff behavior seen in the data. 

Fitting to the dominant van't Hoff behavior ( |55| ) introduces natural scales for energy, m, and 
pressure, C~^, so it is convenient and enlightening to introduce dimensionless variables. With the 
definitions 

/3 = m/3, C{i3) = C0/m)/Co, and F0) = F{l3/m)/mCo, (61) 
the Langmuir curve can be expressed in the dimensionless form 

C0) = mP)e^- (62) 

For consistency with these definitions, the other energy-related functions in the analysis are nondi- 
mensionalized as follows 

G{p) = G{P/m)/m^Co, y = y/m, f{y) = f{my)/Go, ~g{y) = 9i^y)/^Co, (63) 

where f{y) and g{y) are the inverse Laplace transforms of F{(3) and G{(3), respectively. The nat- 
ural scales for energy and pressure also imply natural scales for volume, mGo, and distance, r^, as 
described in the previous section, which motivates the following definitions of the dimensionless cell 
potential versus volume 

X = x/mGo, u{x) = u{mCox)/m (64) 

and radius 

f = r/r^, w{f) = w{rsr) /m. (65) 
Note that x = . With these definitions, the central- well solution takes the simple form, 

u{x) = g~^{x) (66) 

in terms of the dimensionless volume, or 

u,(f) = -1 + ^-i(f3) (67) 

in terms of the dimensionless radius. We now consider various prefactors F{[3) which encode valuable 
information about the energy landscape in various regions of the clathrate cage. 
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6.2 The Interior of the Clathrate Cage 

6.2.1 Power-Law Prefactors 

The simplest possible correction to van't HofF behavior involves a power-law prefactor 

F0) = p-f" for any /x > 0, (68) 

which corresponds to a logarithmic correction on a van't Hoff plot of the Langmuir constant 

logC = /3+(l-^)log(;3) (69) 

as shown in Fig. ^(a). In this case, we have 

f{y) = r-^H{y)/T{n) and ~g{y) = rH{y)/T{^ + 1), (70) 

where T{z) is the gamma function. In general, power-law prefactors at low temperatures signify an 
energy minimum with a simple polynomial shape. 

6.2.2 The Central- Well Solution 

The unique central-well solution is also a simple power law 

n(x) = [r(^+l)x]^/^ (71) 

or equivalently 

u,(f) = -l + r(^ + l)^/^f3/M_ (72) 

The cubic van't Hoff behavior is recovered in the case ^ = 1, as is the (asymptotic) parabolic behavior 
from ( p3D and ( pTj ) in the case /i = 3/2. Because w{f) + 1 oc f^/'^, a power-law correction to van't 
Hoff behavior with a positive exponent (// < 1) corresponds one which is "wider" than a cubic, 
while a negative exponent (/i > 1) corresponds to a potential which is "more narrow" than a cubic, 
as shown in Fig. H(b). On physical grounds, the smooth polynomial behavior described by ( [72[ ) is 
always to be expected near the minimum energy of the cell potential. Therefore, the power-law 
correction to van't Hoff behavior ( |69|) has greatest relevance for low temperature measurements in 
the range /? ^ 1, from which it determines interatomic forces in the interior of the clathrate cage at 
low energies |?i'(r)| ^ 1. 



6.2.3 Non-Central-Well Solutions 



As described above, there are infinitely many non-central-well solutions. One family of solutions of 
the form (p9| ) with A = 1/2 is given by 



w{f) + 1 



[2T{^l + 1) 
[r(^ + i)f3] 



;:31 l/P' 



if < f < 
if f > 2^1^ fn 



(73) 



where = rg/r^ is arbitrary, as shown in Fig. ^(c) for the case = 0.65. These solutions are 
unphysical since they all have cusps at f = 2^/^fo near the outer wall of the cage. However, they 
can still have reasonable behavior near the minimum at fo for certain values of ^, which could have 
experimental relevance for low temperature measurements. Near the minimum, the exact solutions 
(^) have the asymptotic form 



w{f) 



-1 + 



6r(/i + i)f2|f 



1/m 



as r 



(74) 
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which is cusp-Uke for /x > 1/2, but differentiable for < /i < 1/2. For example, the non-central-weh 
has a parabohc shape in the case /i = 1/2, which agrees with the asymptotic analysis in Eqs. (p^)-(|25|) 
when the units are restored, and it has a cubic shape when // = 1/3. On the other hand, in the central- 
well case ^ = 3/2 and fi = 1 correspond to parabolic and cubic minima, respectively. Therefore, 
this example nicely illustrates the difference between the low-energy asymptotics of central and non- 
central-wells described above in section |3|, which would be useful in interpreting any experimental 
Langmuir constant data showing deviations from van't Hoff behavior. 

6.3 The Outer Wall of the Clathrate Cage 

6.3.1 Rational Function Prefactors 

The behavior of the Langmuir curve in the high temperature region /3 = 0(1) is directly linked to 
properties of the outer wall of the clathrate cage, described by the cell potential at high energies w{r)+ 
1 = 0(1). Although this region of the Langmuir curve does not appear to be accessible in experiments 
(see below), in this section we derive exact solutions possessing different kinds of outer walls, whose 
faint signature might someday be observed in experiments at moderate temperatures. In order to 
isolate possible effects of the outer wall, we consider Langmuir curves which are exactly asymptotic 
to the usual van't Hoff behavior at low temperatures with small logarithmic corrections (on a van't 
Hoff plot) at moderate temperatures. These constraints suggest choosing rational functions for F{(3) 
such that F{(3) ~ 1//3 as /3 — > oo. 

6.3.2 Central Wells with Hard Walls 

We begin by considering a "shifted power-law" prefactor 

F0) = 1/0 + a), for any a > (75) 

which corresponds to a shifted logarithmic deviation from van't Hoff behavior, 

logC' = /3-log(l + a//3). (76) 

As shown in Fig. ^(a), this suppresses the Langmuir constant at high temperatures, which intuitively 
should be connected with an enhancement of the strength of the outer wall compared to the cubic 
van't Hoff solution. Taking inverse Laplace transforms we have 

f{y) = e-'^yRiy) and m = (1 - e-"y)H{y)/a, (77) 

and indeed, since g{y) is bounded, all solutions must have a hard wall regardless of whether or not 
the well is central, as described above. For example, the unique central-well solution is 

w{f) = -1 - log(l - ar^) I a for < f < a"^/^ (78) 

which has an outer hard wall at f = as shown in Fig. |6|(b). The solution is also asymptotic 

to the cubic van't Hoff solution at small radii f <C a'^l"^ . Therefore, in the limit a ^ 0, the radius 
of the outer hard wall diverges, and the solution reduces to the cubic shape as the deviation from 
van't Hoff behavior is moved to increasingly large temperatures. Since empirical fitting with Kihara 
potential forms arbitrarily assumes an outer hard wall, this example provides analytical insight into 
the nature of the approximation at moderate to high temperatures, where the Langmuir constant 
should be suppressed according to ([76|). 
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6.3.3 Central Wells with Soft Walls 



Next we consider tlie opposite case of a Langmuir constant whicli is enlianced at high temperatures 
compared to van't Hoff behavior, which intuitively should indicate the presence of a "soft wall", 
rising much less steeply than a cubic function. An convenient choice is 

= - 7^), for any 7 > 0. (79) 

which is analytic except for poles at /3 = ±7 on the real axis. Although this function diverges at 
P = J due to the overly soft outer wall, the corresponding Langmuir curve 

logC' = /3-log[l-(7//3)2] (80) 

shown in Fig. |^(a) could have experimental relevance at moderate temperatures /3 S> 7, if 7 were 
sufficiently small. In this case, we have 

f{y) = cosh{jy)H{y) and g{y) = sinh{-fy)H{y)/j (81) 

which yields the central-well solution 

w{f) = -1 + sinh-^(7f3)/7. (82) 

As shown in Fig. ^(b), this function follows the van't Hoff cubic at small radii f ^ ^~^/^ but "softens" 
to a logarithmic dependence for large radii f ^ 7^^/^. 



7 Interpretation of Experimental Data 

We begin by fitting Langmuir curves, computed from experimental phase equilibria data, an equation 



of state, and reference thermodynamic properties [28| for ethane and cyclopropane clathrate- hydrates 
to the van't Hoff equation 

log C = m/3 + 6 (83) 

using least-squares linear regression. This leads to rather accurate results, as indicated by the small 
uncertainties in the parameters displayed in Table || (63% confidence intervals corresponding to much 
less than one percent error) . The high quality of the regression of log C on /? is further indicated 
by correlation coefficients very close to unity, 0.99650 and 0.99998 for the ethane and cyclopropane 
data, respectively. Using the fitted values for m and Co = e*, the data for the two clathrate- hydrates 
can be combined into a single plot in terms of the dimensionless variables C and /?, as shown in in 
Fig. ^, which further demonstrates the common linear dependence. 

Converting the experimental data to dimensionless variables also reveals that the measurements 
correspond to extremely "low temperatures". This is indicated by large values of /? = m/kT in the 
range of 16 to 24, which imply that kT is less than 6% of the well depth m. As such, physical intuition 
tells us that the experiments can probe the cell potential only very close to its minimum. This 
intuition is firmly supported by the asymptotic analysis above, which (converted to dimensionless 
variables) links the asymptotics of the Langmuir constant for /? ^ 1 to that of the cell potential for 
1^ — fo\ ^ 1- In this light, it is clear that any features of the cell potential other than the local shape 
of its minimum, which are determined by empirical fitting, e.g. using Eq. (^) based on the Kihara 
potential, are simply artifacts of an ad hoc functional form, devoid of any physical significance. 

Since the shape of the potential very close to its minimum should always be well approximated 
by a polynomial (the leading term in its Taylor expansion), the analysis above implies that only 
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simple power-law prefactors to van't Hoff behavior should be considered in fitting low temperature 
data. Therefore, we refit the experimental data, allowing for a logarithmic correction, 

logC = m/3 + 6 + z^log(/?) (84) 

as in Eq. (|69[). The results are shown in Table and the best-fit functions are displayed in dimen- 
sionless form in Fig. |^. In the case of ethane, the best-fit value of ^ = 1 — u corresponds to a roughly 
linear central-well solution w oc r^'^ or a cusp-like non-central-well solution w oc \t — fo\^'^- Although 
these solutions are not physically reasonable, perhaps the qualitative increase in // compared with 
ideal van't Hoff behavior (// = 1) is indicative of a parabolic central well {fi = 3/2). In the case 
of cyclopropane, we have /x = —1.4 it 0.9, which violates the general condition /x > needed for 
the existence of solutions to the inverse problem. If this fit were deemed reliable, then the basic 
postulate of vdWP theory, Eq. (^), would be directly contradicted, with or without the spherical cell 
approximation (see the Appendix B). It is perhaps more likely that the trend of decreasing // < 1 
could indicate a non-central parabolic minimum in the spherically averaged cell potential (/i = 1/2). 

Although it appears there may be systematic deviations from ideal van't Hoff behavior in the 
experimental data for ethane and cyclopropane, 7^ or /i 7^ 1, the results are statistically ambigu- 
ous. For both types of guest molecules, adding the third degree of freedom u substantially degrades 
the accuracy of the two linear parameters m and b, with errors increased by several hundred per- 
cent. Moreover, the uncertainty in z/ is comparable to its best-fit value. Therefore, it seems that 
we cannot trust the results with z/ 7^ 0, and, by the principle of Occam's razor, we are left with the 
more parsimonious two-parameter fit to van't Hoff behavior, which after all is quite good, and its 
associated simple cubic, central-well solution. 

On the other hand, there are different two-parameter fits, motivated by the inversion theory, which 
can describe the experimental data equally well, but which are somewhat more appealing than the 
cubic solution in that they possess a non-vanishing elastic constant (second spatial derivative of the 
energy). For example, the fits can be done using ( |8^ with the parameter f fixed at either 1/2 or 
— 1/2, corresponding to either a non-central or central, parabolic minimum, respectively. The results 
shown in Table |^ reveal that these physically significant changes in the functional form have little 
effect on the van't Hoff parameters m and b = logCo. 

The difficulty with the experimental data as a starting point for inversion is its limited range in 
/? of roughly one decade, which makes it nearly impossible to detect corrections proportional to log /3 
related to different polynomial shapes of the minimum. It would be very useful to extend the range 
of the data, using the analytical predictions to interpret the results. In general, it is notoriously 
difficult to determine power-law prefactors multiplying a dominant exponential dependence, but at 
least the present analysis provides important guidance regarding the appropriate fitting functions, 
which could not be obtained by ad hoc numerical fitting. Moreover, the clear physical meaning of 
the dominant van't Hoff parameters elucidated by the analysis also makes them much more suitable 
to describe experimental data than the artificial Kihara potential parameters. 

8 Summary 

In this article, we have shown that spherically averaged intermolecular potentials can be determined 
analytically from the temperature dependence of Langmuir constants. Starting from the statistical 
theory of van der Waals and Platteeuw, the method has been developed for the case of clathrate- 
hydrates which contain a single type of guest molecule occupying a single type of cage. Finally, 
the method has been applied to experimental data for ethane and cyclopropane clathrate-hydrates. 
Various conclusions of the analysis are summarized below. 

General Theoretical Conclusions 
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• Physically reasonable intermolecular potentials (which are piecewise continuous and bounded 
below) exist only if the Langmuir curve has a dominant exponential (van't Hoff ) dependence 
at low temperatures, lim/3^00 log C//3 = m, with a prefactor F{P) = C{(3)e~"^^ / P which is 
smooth and non-increasing. 

• The slope m of an experimental "van't Hoff plot" of log C versus inverse temperature /5 is 
precisely equal to the well depth, i.e. (minus) the minimum of the potential. This is true 
not only for the spherically averaged cell potential, mm w{r) = —m, but also for the exact 
multi-dimensional potential, min<I>(r) = —m. 

• For any physically reasonable Langmuir curve, the unique central-well potential can be deter- 



mined from Eq. (37) 



There also exist infinitely many non-central- well solutions of the form (|4T| ) , constrained only to 
satisfy Eq. (^). Several classes of such solutions with a central "soft core" (a finite maximum 
at the center of the cage) are described explicitly in Eqs. (^5|)-(|49|). 



• Each one of the multitude of non-central-well solutions with a soft-core typically possesses 
unphysical cusps (slope discontinuities), while the unique central- well solution is a well-behaved 
analytic function. 

• For ideal van't Hoff temperature dependence, C(/3) = Cog"^^, the central- well solution is a 
simple cubic given by Eq. ([5^). The attractive region of the potential has depth m, volume 
mCo, and radius = (3mCo/47r)^/^. Each non-central- well solution for van't Hoff dependence 
has two unphysical cusps, one at the minimum. 

• The experimental signature of a parabolic, non-central-well is a Langmuir curve that behaves 
like C(/3) ~ CqC™^/?^/^ at low temperatures {(3 — > 00), while a parabolic central well corre- 
sponds to C(/3) ~ ae™^/3-i/2. 

• If there is a pure power-law prefactor multiplying van't Hoff behavior C(/3) = Co{mP)^~^e"^^ 
with /U > 0, the central- well solution is also a power-law (^). For certain values of the prefactor 
exponent < fi < 1/2, there are also non-central-well solutions with differentiable minima such 
as (|73|), although such solutions still possess cusps at higher energies. 



• Rational function prefactors multiplying van't Hoff behavior, such as (76) or (|80|), are associated 
with non-cubic behavior at the outer wall of the cage, such as a "hard wall" (|7^) or a "soft 
wall" (|82|), respectively. 



Conclusions for Clathrate-Hydrates 

• Since Langmuir constants must increase with temperature, m > 0, on the basis of general 
thermodynamical arguments, the intermolecular potential must be attractive (with a region of 
negative energy). 

• The depth Wq and radius of the attractive region of the cell potential can be estimated 
directly from experimental data using the simple formulae Wq = —m and = (3mCo/47r)^/^ 
without any numerical fitting. The resulting values for ethane and cyclopropane hydrates are 
consistent with typical estimates obtained by other means. 
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The experimental Langmuir constant data for ethane and cyclopropane clathrate-hydrates is 
very well fit by an ideal van't Hoff dependence, which corresponds to a cubic central well, 

as given by Eq. (|56D. However, the data is also equally consistent with a central parabolic well, 

or various non-central (spherically averaged) parabolic wells. The range of temperatures is 
insufficient to distinguish between these cases. 

Experimental data tends to be taken at very "low" temperatures, kT <C m, which means that 
only the region of the potential very close to the minimum \r — ro\ -C is probed. Therefore, 
only simple polynomial functions are to be expected, and fitting to more complicated functional 
forms, such as the Kihara potential, has little physical significance. 

In practical applications to clathrate-hydrates, the full power of our analysis could be exploited 
by measuring Langmuir hydrate constants over a broader range of temperatures than has 
previously been done. 

The availability of the inversion method obviates the need for empirical fitting procedures p^ , 
p^ ], at least for single-component hydrates in which guest molecules occupy only one type of 
cage. Moreover, the method also allows a systematic analysis of empirical functional forms, 
such as the Kihara potential, which cannot be expected to have much predictive power beyond 
the data sets used in parameter fitting. 

The general method of "exact inversion" developed here could also be applied to other multi- 
phase chemical systems, including guest-molecule adsorption at solid surfaces or in bulk liquids. 
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Appendix A: Inversion of Cohesive Energy Curves for Solids 

The basic idea of obtaining interatomic potentials by "exact inversion" has also recently been pursued 
in solid state physics (albeit based on a very different mathematical formalism having nothing to 
do with statistical mechanics). The inversion approach was pioneered by Carlsson, Gelatt and 



Ehrenreich in 1980 in the case of pair potentials for crystalline metals |14, 15 1. These authors had 
the following insight: Assuming that the total (zero-temperature) cohesive energy E{x) of a crystal 
with nearest neighbor distance x can be expressed as a lattice sum over all pairs of atoms 

XSij j (85) 

where Sij are normalized atomic separation distances, then a unique pair potential (p{r) can be 
derived which exactly reproduces the cohesive energy curve E{x). (Lattice sums also appear in some 
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clathrate- hydrates models [36|, but to our knowledge they have never been used as the basis for an 
inversion procedure.) 

The mathematical theory for the inversion of cohesive energy curves has been developed consid- 
erably in recent years and applied to wide variety of solids |l6|, 0, [l^, 19]. The extension of the 



inversion formalism to semiconductors has required solving a nonlinear generalization of Eq. (|85|) 



representing many-body angle-dependent interactions [18, |T 



^i^) = ^9{xsij)g{xsik)h{dijk) (86) 

ijk 

where F{x) is the many-body energy, h{9ijk) is the energy of the angle between two covalent bonds 
Tij = xsij and fik = xsik, and g{r) is a radial function which sets the range of the interaction. In 
general, (j){r), g{r), and h{9) can be systematically obtained from a set of multiple cohesive energy 
curves for the same material |jl^, The angular interaction can also be obtained directly from 

cohesive energy curves for non-isotopic strains [ST]]. 



Appendix B: Mathematical Theorems 

The first theorem provides necessary conditions on the Langmuir C(/3) so that the cell potential w{r) 
is bounded below and continuous. It also interprets the slope of a van't Hoff plot of the Langmuir 
curve in the low temperature limit as the well depth, under very general conditions. As pointed out 
in the main text, it is convenient to view the inverse temperature /3 as a complex variable. 

Theorem 1 Let w{r) he real and continuous (except at possibly a finite number of discontinuities) 
for r > with a finite minimum, w{r) > Wq = w{ro) > — oo for some ro > 0, and suppose that the 
integral 

POO 

C{l3) = 47r/3 / e-^'^(''V2dr (87) 
Jo 

converges for some (5 = c on the real axis. Then 

C{l3) = /5F(/3)e-"^°^ (88) 

where the complex function F{(3) is 

(i) real, positive and non-increasing on the real axis for (3 > c and 

(ii) analytic in the half plane Re(3 > c. 

If, in addition, the set = {r > 0\wo < w{r) < Wq + e} has nonzero measure for some e = > 0, 
then F{f5) is strictly decreasing on the positive real axis (for [3 > c). Moreover, if has finite, 
nonzero measure for every < e < Eo, then 

lim logC(;3)//3 = -TWo (89) 
where the limit is taken on the real axis. 

Proof: Define a shifted cell potential versus volume, u {^^f^^ = w{r) — Wq. Substituting u{x) for 
w{r) reduces Eq. ( p7[ ) to Eq. (p8|), where 

F{(3) = / e-^"(^)dx. (90) 
Jo 
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Since u{x) > is real, the function F{p) is real and positive for all real (3 for which the integral 
converges. Moreover, for any complex /? and /?' with Re/3 > Re/3' > c, we have the bound 

\F{P)\< e-^^^-"(^)(ix< / e-^^^'-"(^)dx < F(c) < oo (91) 
Jo Jo 



which establishes that the defining integral {p[\ ) converges in the right half plane Re/3 > c and is 
non- increasing on the real axis, thus completing the proof of (i). 

Next let w{r) be larger than its minimum value (but finite), w{ro) < w{r) < oo, on a set Soo 
of nonzero measure, so that < u{x) < oo for the corresponding set of volumes. Then for every 
/3 > /3' > c on the real axis we have 

e-^"(^)d2;< / e-^'"(^)dx. (92) 
On the complement = (0, oo) \ 5*00, either u{x) = or u{x) = oo, which implies 

e-l^^i^)dx= I e-^'"(^)dx. (93) 

CC / CC 

^Prom Eqs. (p^)-(p3D we conclude that F(/3) is strictly decreasing on the real axis. 

Next we establish the low-temperature limit (^). Given < e < eo, we have the following lower 
bound for any /3 > c on the real axis: 

e'PF{l3) = /°°e-^["(^)-^](ix > e'f^'^ [ dx + [ e-'^^^'^'^^-'Ux > e^f^l'^ [ dx. (94) 
Jo Js,/2 Jsif^ Js,/2 

Combining this with the upper bound, F{(3) < F{c) < oo, we obtain 

e-'/'^M, < F{I3) < F{c) (95) 

where M^: = Jg ^ dx is a finite, nonzero constant (because is assumed to have finite, nonzero 
measure). Substituting Eq. ( |88[ ) in Eq. (95), we arrive at 

log/3 + logM, -e/3/2-/3'«;o < logC7(/3) < log^ + logF(c) - ^u;^ (96) 

which yields 

-Wo-e/2 < linilogC(/3)//3 < -Wo- (97) 

The desired result is obtained in the limit e — > 0. 

Finally, we establish the analyticity of F{P) in the open half plane Re/3 > c by showing that its 
derivative exists and is given explicitly by 

F'{0) = - / e-^"(^')n(x)dx. (98) 



This requires justifying the passing a derivative inside the integral (20), which we have just shown to 
converge for Re/3 > c. Using a classical theorem of analysis [^], it suffices to show that the integral 
in ( |98[ ) converges uniformly for Re/3 > c + e for every e > because the integrand is a continuous 
function of (3 and x. (The possibility of a finite number of discontinuities in u{x) is easily handled by 
expressing ( p8| ) as finite sum of integrals with continuous integrands.) It is a simple calculus exercise 
to show that te^* <l/e, and hence 

< ^ (99) 



ee 
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for all real t > 0. This allows us to derive a bound on the "tail" of the integral ([Qq): 

roo . roo 1 roc 



e 

X ' Jx Jx ' ee Jx 

which is independent of /?. This uniform bound vanishes in the hmit X — > oo because it is proportional 
to the tail of the convergent integral defining -F(/3), which completes the proof. □ 

The proof of Theorem 1 does not depend in any way on the dimensionality of the integral and 
thus can be trivially extended to the general multi-dimensional case of vdWP theory without the 
spherical cell approximation. 

Theorem 2 Let <I>(r, 6, cj), a, ^, 7) > $(ro, 9o, (po, cto, Co, 7o) = Wo be real and continuous, and suppose 
that the integral 

C{(3) = ^ / e-^^^~'''^^^^'''^'^\'^sinesinidrded(^dadid-^ (101) 



87r2 jv 

converges for some (3 = c (real). Then all the conclusions of Theorem 1 hold. 

The six-dimensional integral ( |101| ) of Theorem 2 does not present a well-posed inverse problem for 
the intermolecular potential However, the spherically averaged integral equation ( ^7|) of Theorem 
1 can be solved for the cell potential w{r) for a broad class of Langmuir curves C(/3) specified in the 
following theorem. The proof is spread throughout section ^ of the main text. 

Theorem 3 If the inverse Laplace transform f{y) of F{j3) exists and is nondecreasing and non- 
constant for y > 0, then there exist a unique central-well solution (ro = Oj and infinitely many 



non- central-well solutions (0 < ro < 00) to the inverse problem (11). If f{y) is also continuous, then 



the central-well solution is the only continuously differentiable solution. 

Finally, we state sufficient assumptions on F{I3) to guarantee the assumed properties of f{y). 
In light of the necessary condition that F{P) be analytic the right half plane Re/3 > c, the defining 
contour integral for /(y) 

1 rc'+ioo 

f{y) = — / e^^F(/3)d/3 (102) 



C —too 
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must converge for any c' > c. By closing the contour in the left half plane, it can be shown that a 
sufficient (but not necessary) condition to ensure the assumed properties of f{y) is that F{j3) decay 
in the left half plane (limp^oo W{p^^^) \ = for tt/2 < 9 < 37r/2) and have isolated singularities only 
on the negative real axis or at the origin with positive real residues. The particular examples of F{P) 
considered in section ^ satisfy these conditions, but the weaker assumptions above regarding f{y) 
suffice for the general derivation in section ^. 
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Figure 1: Exponential fits of Langmuir constants over the measured tem- 
perature range plotted with linear axes for (a)-(b) cyclopropane and (c) 
ethane clathr ate- hydrates. An enlargement of the high temperature data for 
cyclopropane is shown in (a). The experimental data is taken from Ref. p8| . 
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Figure 2: Exponential dependence with inverse temperature of experimental 
Langmuir curves from Fig. || plotted with log-linear axes for (a) ethane and 
(b) cyclopropane clathrate hydrates. Straight lines indicate pure van't Hoff 
behavior. 
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Figure 3: (a) Sketch of a central-well solution, where x is the scaled volume 
of interaction and u{x) is the spherically averaged cell potential with inverse 
v{y) = u~^{y). (b) Sketch of a non-central- well solution composed of a non- 
increasing function u^(x) and a nondecreasing function u^{x) joined at a 
minimum of zero at Xq, along with a possible hard core at xi and hard wall 
at X2- The two branches V-{y) and v+{y) of the multi- valued inverse cell 
potential v{y) are also shown, along with other variables defined in the text. 
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Figure 4: Analytical cell potentials for the ethane clathrate-hydrate which 
exactly reproduce the experimental data in Fig. 0(a). The unique central- 
well solution (56) is indicated by a solid line, while a family of non-central- 
well solutions with soft cores ( |59| ) is also shown as dashed lines with cusp- 
like minima at = 0.2, 0.3, 0.4 A. Each of these solutions also has a cusp 
at r = 2^/'^ro, where the energy is the same as the central maximum, and 
beyond this distance joins the central-well solution. 
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Figure 5: Exact inversion of Langmuir curves with power-law corrections to 
van't Hoff behavior in terms of dimensionless variables, as in Eq. (|69|). (a) 
Plots of C = CjCo versus /3 = mjkT for the cases /i = 1/2, 1,3/2. (b) The 
corresponding (unique) central- well potentials plotted as w = w/m versus 
f = r/(3mCo/47r)^/'^. (c) Examples of soft-core non-central- well solutions 
of the form ([7^) with an arbitrarily chosen minimum at f = 0.65, which all 
have cusps at f = 2^/3 (0.65) « 0.819. 
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Figure 6: Exact inversion of Langmuir curves of the form C = e^/(l + 
a/ (3 — (7//?)^) in terms of the dimensionless variables defined in Fig. ^. (a) 
Langmuir curves in this class of functions have anomalous high temperature 
(small (3) behavior but are asymptotic to the van't Hoff curve (a = 7 = 0). 
(b) The corresponding central-well solutions depart from the cubic van't 
Hoff curve at large radius and energy, indicating different properties at the 
boundary of the clathrate cage: "hard walls", if a > and 7 = 0, or "soft 
walls", if a = and 7 > 0. 
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Figure 7: Collapse of the experimental Langmuir curves for ethane and 
cyclopropane hydrates from Fig. ^ plotted in terms of the dimensionless 
variables C = C/Cq versus (3 = m/kT, where Co and m are obtained by 
fitting each data set to logC = m/kT + log Co- Ideal van't Hoff behavior 
C = exp(/3) is shown as a solid line. Fits including power law corrections as 
in Eq. ( |8^ are also shown as the dotted and dashed lines (which are very 
close to the van't Hoff line). 
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Table 1: Linear regressions of the experimental Langmuir constant data 
for ethane and cyclopropane clathrate-hydrates on the form log C = m(3 + 
b + z^log(/3), where b = log Co. Errors reflect 63% confidence intervals. The 
rows where 1^ = 0, 1/2, or —1/2 correspond to two-parameter fits with held 
constant. 



Guest Molecule 


m (kcal/mol) 


log Co (atm ^) 


V 


Ethane 


9.422 ± 0.043 


-14.561 ± 0.081 







9.180 ± 0.044 


-14.419 ± 0.082 


1/2 




9.664 ± 0.043 


-14.703 ± 0.080 


-1/2 




10.52 ± 0.85 


-15.2 ± 0.50 


-2.3 ± 1.8 


Cyclopropane 


10.594 ± 0.012 


-15.474 ± 0.022 







10.335 ± 0.011 


-15.302 ± 0.021 


1/2 




10.853 ± 0.012 


-15.646 ± 0.024 


-1/2 




9.36 ± 0.47 


-14.66 ± 0.31 


2.37 ± 0.90 
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